Phase separation in peptide aggregation processes - 
Multicanonical study of a mesoscopic model 
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We have performed multicanonical computer simulations of a small system of short protein- 
like heteropolymers and found that their aggregation transition possesses similarities to first- 
order phase separation processes. Not being a phase transition in the thermodynamic sense, the 
observed folding-binding behavior exhibits fascinating features leading to the conclusion that 
the temperature is no suitable control parameter in the transition region. More formally, for 
such small systems the microcanonical inteipretation is more favorable than the typically used 
canonical picture. 

1 Introduction 

Folding-binding and docking processes between proteins are significant for catalysis, trans- 
port, and cell stabilization in biological systems. Also, gene replication and expression are 
impossible without defined binding mechanisms of molecules. However, the mutual influ- 
ence of proteins on each other can also result in refolding of proteins (which often leads to 
the loss of their functionality and thus biological activity) or cluster formation. In the latter 
case, proteins self-assemble and form aggregates. The effects of plaque can be disastrous 
and cause heavy diseases: First, the assembled proteins lose their individual functionality 
and second, in the passive case, the aggregates hinder transport and signal exchange pro- 
cesses which are significant for the life of cells. In an active process, specific aggregates 
might be able to bind to cell membranes and to change the membrane morphology, e.g., by 
forming pores. In the amyloid hypothesis for the onset of Alzheimer's disease, for exam- 
ple, aggregates of A/3 proteins are believed to form pores in membranes of neuron cells, 
thus opening ion channels for neurotoxic calcium. ' 

We focus here on thermodynamic properties of the aggregation transition of small pep- 
tides. For this purpose, a simple hydrophobic -polar aggregation model is introduced and 
employed in a multicanonical study of a few short heteropolymers.^' '' 

2 Aggregation model 

For the aggregation study, we extend the AB model^ by an additional interchain interaction 
between the M heteropolymers. As in the single-chain model, which has proven to be quite 



useful in qualitative studies of tertiary folding behavior, only two types of amino acids are 
considered: hydrophobic residues (A) which avoid contact with the polar environment and 
polar residues (B) being favorably attracted by the solvent. The single-chain energy of the 
/ith heteropolymer (/i = 1, • ■ • , M) composed of monomers is given by^ 



where < "di ^ < it denotes the virtual bending angle between the monomers i^, + 1, 
and + 2. Not discriminating nonbonded interactions between monomers of the same or 
different polymers, our aggregation model reads^ 



where /i, label the M polymers interacting with each other, and ifj,,ji, index the 
monomers of the respective /ith and z^th polymer. The nonbonded interresidue pair po- 
tential *(r.vj„;cr.v,crjj = ^r~f^ - Ciai^,aj^)r~]j depends on distance r,^j„ and 
residue type <Ti ^ — A, B. The long-range behavior is attractive for like pairs of residues 
[C{A,A) = 1,"C(B,S) = 0.5] and repulsive otherwise [C{A,B) = C{B,A) = -0.5]. 
The length of all virtual peptide bonds is unity. In this short note, we focus on a sys- 
tem of two identical chains with the Fibonacci sequence AB2AB2ABAB2AB, where the 
single-chain properties are known. ^ Our primary interest is devoted to the phase behavior 
of the system and for this purpose, the density of states g{E) is a suitable quantity that we 
obtained by means of multicanonical computer simulations.^ 

3 IVIicrocanoiiical vs. canonical view 

The Hertz definition of the entropy is given by S{E) = ks In r(£'), where ks is the 
Boltzmann constant (ks = 1 in our simulations) and r(£') = /| dE' g{E') (where 
i?min is the ground-state energy) is the phase space volume. In Fig.[Tl S{E) is shown for 
our two-peptide system. Interestingly, in the energy region between -Eagg and -Eftag the 
entropy exhibits a convex behavior, which is a strong indication for surface effects within 
this small system.** Also shown in Fig.[T]is the corresponding concave hull Hs{E), i.e., the 
Gibbs construction. The surface entropy, defined as AS{E) — Ti,s{E) — S{E) is maximal 
at the energy E^cp- The reason for the nonvanishing surface entropy is that the transition 
between the fragmented, i.e., separated chains, and the formation of a joint aggregate is 
a process with phase separation which is "delayed" due to steric surface effects reducing 
the entropy of the total system. Since entropy reduction is only achieved by additional 
energy consumption, the surprising side effect is that in the transition regime the aggregate 
becomes colder with increasing system energy. This is verified by considering the caloric 
temperature which is defined via T~^{E) = dS{E) / dE, also shown in Fig.[T] Actually, 
in the transition region, T^-^{E) bends back with increasing energy. 

Consequently, there is no unique mapping between temperature and energy in the tran- 
sition region (or more precisely, within the bounds T^^ and T> ^ indicated in Fig.!!). Thus 
the temperature should not be considered as a suitable external control parameter. From 
a statistical point of view this means that for transitions with phase separation in small 
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Figure 1. Microcanonical Hertz entropy S{E), concave Gibbs hull HsiE), and inverse caloric temperature 
T~^{E) as fimctions of energy. The phase separation regime is bounded by Eagg and iJfrag; the tem- 
perature region, where temperature is no suitable external control parameter and the canonical interpretation 
breaks down, ranges from TJ^ to . The slope of the Gibbs huU defines the aggregation temperature, 
Ta-^g = d-Hs{E)/dE = const. 

systems a microcanonical interpretation is preferred over the typically used canonical for- 
malism. Since the backbending effect in the peptide aggregation process is a real physical 
effect, it should also be accessible to experimental verification, as it has indeed already 
been observed, for example, in experiments of sodium cluster formation processes.^ 
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